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Theory of heavy quark energy loss 
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We briefly review some of the models and theoretical schemes established to describe 
heavy quark quenching in ultrarelativistic heavy ions collisions. Some lessons are derived 
from RHIC and early LHC data, especially as for the contraints they impose on those models. 

<N ■ 

§1. Introduction 

o 

Jet-quenching of light hadrons is one of the key features observed in ultra- 
relativistic heavy ions collisions (URHIC) performed at RHIC (y^s = 200 AGeV). 
It has been recently confirmed in yfs = 2.76 j4GeV lead-lead collisions at the LHC 
(see^ 1 and references therein). This observation is usually interpreted as the conse- 
quence of the energy loss which affects the leading parton during its passage through 
a hot deconfined medium and is therefore crucial for scrutinizing its properties. In 
recent years, several theoretical schemes, based on the eikonal approximation, have 
^p -1 , been developed in order to describe this energy lossP'®'®'^ Numerical evaluations 

C^ij of these approaches turned out to be able to reproduce the rather flat momentum 

dependence of the nuclear modification factor (Raa) observed at RHIC for pions 
at large px- The driving parameters needed to obtain quantitative agreement P 
however, overshoot pQCD predictions, sometimes by a factor IcPI This might be 
J> ■ taken as an indication that the fundamental theory is not able yet to describe the 

experimental results without introducing ad hoc parameters. 

This is one of the motivations for addressing on the same footing the quenching 
■^J- ■ of jets consisting of leading heavy quarks (HQ), as they might help to better con- 

strain the models. In this respect, the guiding concept is the so-called mass-hierarchy 
AE(g) > AE(q) > AE(c) > AE{b). The first inequality stems from the respective 
SU(3) Casimirs of the gluons and quarks. The second is generally attributed to 
the dead cone effeclP in the context of radiative energy loss but is present when- 
ever the higher mass of the parton implies a reduction of the formation time and 
^ ■ hence of the radiated field. It is also found in collisional energy lossP Before the 

advent of the HQ data at RHIC, it was even advocated^ 1 that HQ jets might be un- 
quenched, but early Raa and elliptic flow (V2) data of non-photonic single-electrons 
(NPSE) revealed that also HQ were strongly quenched in those collisions, nearly as 
much as light onesP^ Since then, some schemes designed for light quarks have been 
extended to heavy quarks^ 1 JEDE) with the main conclusion that a common quan- 
titative agreement between ^^(71") and i?yiA(NPSE) can only be achieved if those 
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*' This large discrepancy should be taken with a grain of salt, as the optimal parameter vastly 
depends on the way the underlying medium is described in the model (seeP for a recent discussion, 
as well as the contribution of S. Bass to these proceedings). 
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NPSE stem exclusively from c-quarks. This is, however, incompatible with recent 
STAR measurements in p — jflS and FONLL calculations® which indicate a signif- 
icant component of 6-quarks for pt > Pt,ctoss ~ 5 GeV/c. Faced with this puzzle 
of lack of coupling of HQ with the QGP, several dedicated models have nourished 
in the literature. In this contribution, we will review them shorthF^l Then we will 
provide a tentative explanation why several models with various physical inputs are 
able to cope with data. Finally, we address the benefits of LHC and conclude. 

§2. Models at RHIC: fragility and robustness 

2.1. The two faces of heavy quark energy loss 

The first approach of HQ propagation in a hot medium was based on the Fokker- 
Planck equation with drag and diffusion coefficients evaluated from collisional energy 
loss only.® This type of approach has survived up to now,® although the physical 
content of the basic interaction has been extended by several authors (presented by 
increasing order of "strong coupling" content): A. Peshier® underlined the role of 
genuine running a s in collisional energy loss. Calculations implementing such run- 
ning feature®® 1 indeed find a smaller discrepancy with the iiU^NPSE), sometimes 
at the price of an additional cranking factor of the order of 2 — 4. Van Hees and Rapp 
were the first who have investigated the effect of possible heavy-light quark bound 
states in the s-channel.® Later, they have proposed a T-matrix description of HQ 
diffusion® based on a bona-fide generalization of the static Q — Q potential evalu- 
ated through lattice calculations. In both cases they obtain a good agreement with 
Raa as well as with t>2 data. Akamatsu et al.® have implemented the Langevin evo- 
lution of HQ in a hydrodynamical medium resorting to drag and diffusion coefficient 
evaluated in the strong coupling limit through AdS/CFT correspondence. Strictly 
speaking, all these FP/Langevin treatments applies for pt < m. In this regime, HQ 
indeed behave as heavy particles surrounded by light degrees of freedom and their 
thermalization time oc m differ significantly from that of light quarks. 

For ultra-relativistic HQ (pt 3> m), the radiative energy loss becomes the dom- 
inant mechanism. The results of most of the eikonal schemes mentioned above is 
the probability of radiation which presents significant fluctuations. In principle, this 
feature invalidates the FP treatment of HQ propagation at high energy. In this 
regime, the mass of the quark acts mostly as a collinear regulator. For the most en- 
ergetic case (E 3> m), in-medium formation time of the high energy gluons (^Ju/q) 
exceeds the path length L. Then this scale regulates the radiation spectrum. As 
a consequence the average energy loss AE oc L 2 and the mass merely appears (if 
at all) through a logarithmic factor.®® 1 it should be noted that the conversion of 
heavy quarks into heavy mesons differ as well in these two regimes: for px < m, 
coalescence dominates because the probability to pick up a light quark from the 
medium at small relative velocity is large. For px S> m, this probability is close to 
and fragmentation becomes the dominant mechanism, with, however non trivial 
consequences on the Raa due to the presence of possible in-medium bound states.® 



See f.i.^J for an extended review 
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Fig. 1. Raa of non photonic single electrons 
for central Au-Au collisions at RHIC en- 
ergy; see^ for details. 



2.2. Model fragility 

Although the approaches presented above differ vastly w.r.t their physical as- 
sumptions they are all able to cope with the NPSE Raa data measured by RHIC ex- 
periments, at the price of a rescaling of the coupling parameter. 
This questions our ability to achieve ro- 
bust conclusions w.r.t. the basic mech- 
anism at hand and hence the QGP 
properties. IrJ23 we have implemented 
collisional as well as radiative energy 
loss in the same numerical framework 
(MCa s HQ) and confirmed this "model 
fragility" (see fig. [TJ. A finer analy- 
sis reveals that the Raa observable is 
mostly sensitive to the energy loss spec- 
trum P{oo) at low values of the energy 
loss uj. As pointed out by Baier et alp^ 1 
this is due to the rather stiff initial pr- 
distribution which makes the sequence 

of many small losses more probable than a single process involving large energy loss. 
This explains why both spectra (see^J for an illustration) lead to similar quenching 
although the average AE ra _d > AE co \. 

2.3. "Robust" contact with lattice calculations 

According to our analysis, all HQ observables at RHIC can be explained with a 
simple AE oc L law, an observation in favor of local processes. 
On fig. [21 we show the drag coeffi- 
cient A from the different microscopic 
energy-loss models implemented in our 
MCa^HQ framework. Although these 
models vary largely, their values for A 
nicely approach a rather unique value 
f° r Pt ^ 5 GeV/c once the coupling 
parameter (here the interaction rate) is 
rescaled in order to match RHIC NPSE 
data. From this we extract "robust" 
values of the relaxation coefficient rj = 

linip^o of °- 45 ± °- 15 c / fm at T = 
200 MeV and 0.75 ± 0.25 c/fm at T = 
300 MeV. This yields a spacial diffusion 
coefficient D s of 2itTD s = 1.9 ± 0.5 at 

T = 200 MeV and of 2ttTD s = 2.55 ± 0.65 at T = 300 MeV, for both c and b 
quarks, in quite good agreement with recent lattice calculations!^ This first suc- 
cessful contact with lattice should encourage us to seek for alternative observables 
able to discriminate between various models for px > 10 GeV/c. 
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Fig. 2. Drag coefficient of various models in- 
corporated in MCa s HQ (running a s colli- 
sional, incoherent and coherent radiation), 
with coupling parameter rescaled in order 
to match the NPSE RHIC data. 
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2.4. The role of the elliptic flow 

In principle, the vi observable helps in constraining the models. At low and 
intermediate pr, it reflects the collectivity acquired by heavy quarks and develops 
constantly with time up to the end of the transition! ® n [ s thus more sensitive 
to the QGP evolution than the Raa which saturates earlier. Nevertheless, it has 
been shown recently® that two energy loss models with drag factor differing by a 
factor as large as two were both able to reproduce the experimental Raa and v<i 
once they are imbedded in different QGP evolution models chosen by the respective 
authors.®'® This clearly points towards the need of performing joint analysis of 
bulk QGP properties and HQ observables to achieve significant progress, as initiated 
in® for the case of light hadrons. At large px, one expects some V2 as well, due to 
the path length difference along both principal directions of the QGP, understood as 
a source of quenching.® Although this observable could be useful to assess the path 
length dependence of HQ quenching (see® for prediction within the ASW model), 
present px range available at RHIC seems to offer little discriminating power for this 
purpose. 



§3. Benefits from early LHC 



URHIC performed at LHC offer larger possibilities to discriminate between 

various models and theoretical schemes, due to wider pt range and due to the 

possibility to measure D and .B-mesons observables separately from the very first 

runs on. Theoretical predictions®-® for the R A a of D -mesons range between 0.2 

and 0.4 for pr > 10 GeV/c. First ALICE results on D mesons in Pb-Pb colli- 
M 



sions^a confirmed the HQ quenching found at RHIC, with Raa ~ 0.3 ± 0.15 for 
p T G [5 GeV/c, 10 GeV/c]. 
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Fig. 3. Quenching of D mesons in PbPb collisions at 2.76 TeV: preliminary ALICE datcP* 
collisional + radiative energy loss implemented in MCa s HQ. 



On fig. [3l we provide the comparison between data and calculations performed 
with the same elementary reaction cross section as for the RHIC energy. We only 
adapted the plasma expansion (simulated with the hydro dynamical code of Kolb- 
Hein^SJ) to obtain = 1600 at the chemical freeze out and the initial px- 

distribution (taken according to FONLL 1.3.2), which is harder than at the RHIC 
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energj0 A satisfactory agreement is obtained although the p^-range of the present 
data is not sufficient to disambiguate between collisional and radiative energy loss. 
Apart from the px dependence of Raa and V2, the various models exhibit a rather 
rich mass-dependence which can be probed by addressing simultaneously the light 
vs heavy or c vs b observables. The quenching of B mesons has not been mea- 
sured directly yet but can be bona fide assimilated to the quenching of non-prompt 
J/ip measured by CMSj^ 1 this is legitimate due to the flat shape of Raa at high 
Pt- On fig. HI we show the comparison between these data and some collected 
predictions P^'^ 1 Most of the predictions show a lack of quenching as compared to 
the data. If confirmed, this would indicate that history might repeat itself (large 
quenching of c quark unpredicted and observed at RHIC; large quenching of b quark 
unpredicted and observed at LHC) and would question our present "understanding" 
of HQ propagation in hot media. On fig. [5] (right), we present various predictions for 
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Fig. 4. Left: quenching of B mesons in PbPb collisions at 2.76 TeV (identified to the quenching of 
non-prompt J/ip measured by CM£p^ (preliminary data) vs various model predictions. Right: 
Rcb predictions vs experimental ratio extracted from ALICE data on D-mesons and CMS data 
on non-prompt J /if) in central PbPb collisions (black disk).; caution: the prediction of Greco et 
al. was made for minimum bias collisions. 



the Rcb ratio (defined as Raa(D) / Raa(B)) as a function of pr, together with the 
experimental point obtained by Raa{B) = Raa (non — prompt J/tp) measured by 
CMS. Although the rather large experimental errors prevent us presently from any 
firm conclusion, models incorporating rather large mass dependence like AdS/CFT 
(drag coefficient oc m _1 ) seem to be disfavored by the data, although the analog of 
the detailed in-medium evolution oP> is lacking at LHC. 



§4. Conclusions and prospects 



We have argued that the limited dynamical range of observables at RHIC hin- 
ders the discrimination between the various models proposed for HQ evolution in 
hot media. We have shown that it was nevertheless possible to extract some basic 
properties of the QGP-HQ interaction - like the drag coefficient at low momentum 

*' This explains why the medium appears as less opaque to HQ propagation although it is 
denser. 
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- in a rather "robust" way, once those models are rescaled in order to match Raa 
and V2 data. These values are in good agreement with present lattice calculations. 
We have discussed the recent LHC results pertaining to .D-mesons and non-prompt 
J/tjj from B meson decay. In particular, we found that most of the predictions seem 
to lack quenching for the 6-quark. Whether it is just a question of fine-tuning of the 
present models or the sign of a more fundamental misunderstanding of the physics 
requires refined data and access to the full pt range of observables like Rcb- They 
will come in the next few years. Finding Raa(q) ~ Raa(c) ~ Raa{o) at mid pt 
could be the sign that the gluon formation-time which is found to be an increasing 
function of the radiator's inverse mass^ZJ could be bounded from above by another 
scale of the problem as for instance the gluon absorption length in the mediumPS 
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